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Abstract —This paper considers a downlink non-orthogonal 
multiple access (NOMA) system where the sonrce intends to 
transmit independent information to the users at targeted data 
rates under statistical channel state information at the transmit¬ 
ter. The problem of outage balancing among the users is studied 
with the issues of power allocation, decoding order selection, and 
user grouping being taken into account. Specifically, with regard 
to the max-min fairness criterion, we derive the optimal power 
allocation in closed-form and prove the corresponding optimal 
decoding order for the elementary downlink NOMA system. By 
assigning a weighting factor for each user, the analytical results 
can be used to evaluate the outage performance of the downlink 
NOMA system under various fairness constraints. Further, we 
investigate the case with user grouping, in which each user group 
can be treated as an elementary downlink NOMA system. The 
associated problems of power and resource allocation among 
different user groups are solved. The implementation complexity 
issue of NOMA is also considered with focus on that caused 
by successive interference cancellation and user grouping. The 
complexity and performance tradeoff is analyzed by simulations, 
which provides fruitful insights for the practical application of 
NOMA. The simulation results substantiate our analysis and 
show considerable performance gain of NOMA when compared 
with orthogonal multiple access. 

Index Terms —non-orthogonal multiple access, statistical chan¬ 
nel state information, outage probability, fairness, optimization. 


I. Introduction 

HE concept of superposition coding (SC) is originally 
proposed for broadcast channels (BCs), where the dis¬ 
parity in the channel qualities of the users due to the near- 
far effect and the random fading nature of the radio channels 
can be exploited as a new degree of freedom for potential 
performance gains IT]. Recently, SC has received renewed 
attentions for its large potential in throughput and user fairness 
enhancement 12], IS]. In particular, non-orthogonal multiple 
access (NOMA) by using SC at the transmitter and successive 
interference cancellation (SIC) at the receiver has been widely 
studied and is recognized as a promising candidate for future 
5G systems. In l4|, Q, the problem of multi-user scheduling 
in NOMA is investigated. In fb], I?] and references therein, 
system-level performance evaluation of NOMA is conducted 
from various perspectives. All these works show considerable 
performance gain of NOMA when compared with orthogonal 
multiple access (OMA). 
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In NOMA systems, independent signals dedicated to dif¬ 
ferent users are transmitted concurrently in the same time- 
frequency channel, which causes strong inter-user interference 
when decoding at the receiver sides. To attain the potential 
performance gain of NOMA, advanced reception technique 
that can distinguish the overlapped signals is required la. SIC 
is the mostly considered for its low complexity and simplicity 
of application in existing systems El, 13. However, SIC gives 
rise to the problem of decoding order selection at the receivers, 
which complicates the performance optimization of NOMA, 
since the decoding order selection problem is of combinatorial 
nature and difficult to solve in many cases of interest. Most of 
the existing work on NOMA assume that instantaneous chan¬ 
nel state information (CSI) is available at the transmitter a- 
0. In this case, the fading BC (which models the downlink 
NOMA system from an information-theoretic viewpoint) can 
be treated as multiple parallel degraded BCs m-in]; hence, 
only the decoding order selection problem of the degraded 
BC needs to be considered. The capacity-achieving decoding 
order of a degraded BC is already very clear in the literature 
HI. However, when CSI is unavailable at the transmitter, the 
problem becomes rather complicated. In this case, only the 
outage capacity region of the fading BC is solved 03, ESI. 
Interested readers can refer to lfT6ll - IT3 and references therein 
for the information-theoretic work on the fading BC without 
transmit CSI. 

The lack of transmit CSI is a non-trivial case of practical 
interest. This may happen when the feedback links are limited 
or for other reasons. In contrast, statistical CSI is easier to 
obtain, especially in NOMA systems where multiple users are 
involved in a transmission. Hence it is improtant to investigate 
the performance of NOMA under statistical CSI. We focus on 
the downlink scenario where the source transmits to each user 
at a targeted data rate, for which the outage probability is an 
appropriate performance metric ll20ll . EH- In ED, the outage 
performance of the downlink NOMA system was analyzed for 
preassigned power allocation and decoding order while not the 
optimal one. The work in ll23l investigated the impact of power 
allocation on the fairness outage performance of the downlink 
NOMA system. However, they did not consider the decoding 
order selection problem, only the optimal power allocation for 
a preassigned decoding order was solved. 

In this paper, power allocation and decoding order selection 
are jointly studied for the downlink NOMA system to balance 
the outage probabilities of the users. As in E3 . the max-min 
fairness criterion is considered. Moveover, to facilitate more 
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flexible fairness modeling, weighting factors are assigned for 
the users to indicate the relative importance of their success 
probabilities. The purpose of this work is two folds, one is 
to provide a theoretical tool for evaluation of the fairness 
outage performance of NOMA in the downlink scenario, the 
other is to analyze the complexity and performance tradeoff 
of NOMA under the investigated system setup. The fairness 
enhanced nature of NOMA has been extensively studied ll24l . 
Ea. However, these work concentrated on the achievable rate 
performance with the assumption that CSI is available at the 
transmitter, fairness with regard to the outage performance is 
not considered. The second purpose is motivated by the fact 
that the complexity of SIC scales at least linearly with the 
number of the users m, m, M, which forms the main 
obstacle to the practical application of NOMA. We investigate 
the effectiveness of user grouping in combating this problem. 
Specifically, we can divide the users into multiple groups and 
schedule only one group of the users at each time. Obviously, 
the complexity of SIC depends on the size of the groups. Also, 
the power and resource allocation among different user groups 
as well as the user grouping algorithm will introduce new 
complexity issues. All these aspects are taken into account 
when analyzing the complexity and performance tradeoff in 
our work. The main contributions of this paper are as follows: 

• Solving the joint power allocation and decoding order 
selection problem analytically for outage balancing in the 
elementary downlink NOMA system where all users are 
scheduled simultaneously (i.e., without user grouping). 
Note that in ll2^ . the power allocation problem was 
investigated for a given decoding order only; moveover, 
only numerical solution of the optimal power allocation 
was obtained by using an iterative search algorithm which 
introduces nontrivial computational complexity. 

• Investigating user grouping in downlink NOMA as well 
as the corresponding inter-group power and resource allo¬ 
cation problems. A low-complexity algorithm is provided 
to obtain the optimal resource allocation among the user 
groups to balance the outage probabilities of the users. 

• Conducting and comparing the simulations with different 
types of resource allocation (continuous, discrete, or 
without) and user grouping (random, optimal, or without) 
to investigate the complexity and performance tradeoff of 
the downlink NOMA system. It is demonstrated that user 
grouping is an effective method in reducing the complex¬ 
ity of NOMA, while causes only moderate performance 
degradation. 

The remainder of this paper is organized as follows. Section 
In] introduces the system model and the optimization problems. 
In section Uni the power allocation and decoding order selec¬ 
tion problems are investigated for the elementary downlink 
NOMA system. The case with user grouping is considered in 
Section|IVl Simulation results are given in SectionlV] Also, the 
complexity and performance tradeoff of NOMA is discussed. 
Section IVTl concludes this paper. 

Notations: Throughout this paper, ]E( 7 ) denotes the expec¬ 
tation of the random variable 7 and we use 7 |a=ao denote 
the value of the variable 7 when a = oq. We denote by Pr(0) 


the probability of the event O and by O the complementary 
event of O. 

II. System Setup and Problem Formulation 

Consider a downlink NOMA system with a source node and 
K destination nodes or users, Uk, fcG/C = {l,2,-' - ,Ar}. All 
the nodes are equipped with a single antenna. The source has 
the mission of delivering mutually independent information to 
the users at a targeted data rate, denoted by for each user 
Uk, k £ 1C. Unlike in conventional OMA systems where each 
channel block can be used by at most one user, in NOMA more 
than one users can be scheduled simultaneously in the same 
channel block. Specifically, the independent signals dedicated 
to the scheduled users are combined at the source using SC and 
then transmitted to the users. At the user sides, SIC is adopted 
to extract the desired information from the combined signal. It 
is assumed that the channels between the source and the users 
undergo independent Rayleigh fadings that are constant over 
one channel block while vary independently from block to 
block (i.e., block fading). Moveover, as in ||23]| . it is assumed 
that perfect CSI is known to the appropriate receivers while 
only statistical CSI is known to the source. 

We consider both the cases with and without user grouping. 
For the case without user grouping, all the users are scheduled 
simultaneously in each channel block, which is termed the 
elementary downlink NOMA system in this paper. In this case, 
the signal received by user Uk can be written as 

Uk = s/lhhkX + Zk,k £ K., (1) 

where x is the signal transmitted by the source, hk is the 
normalized Rayleigh fading coefficient between the source and 
user Uk, Hk is the average channel gain from the source to 
user Uk, and Zk is the additive white gaussian noise with zero 
mean and variance Nq at user Uk- It should be noted that the 
transmit signal x is a superposition of K independent user- 
dedicated signals, which goes 

X = ^ s/akPxk, (2) 

kelC 

where the signal Xk contains the information required by Uk 
and satisfies E(|xfcp) = 1 , P is the short-term transmit power 
constraint, and ak is the power allocation factor (PAF) for 
Uk denoting the proportion of the transmit power allocated to 
Xk- For notation, use 7 ^ = PHk\hk\'^/N q and F^ = ^( 7 ^) 
to denote the channel SNR of user Uk and its mean value, 
respectively. Thus, jk follows an exponential distribution with 
parameter l/F^. 

When decoding at each user Uk, k £ 1C, the desired signal 
Xk is interfered by the other users’ signals. SIC will be 
carried out at the users to mitigate the negative effect of the 
inter-user interference. The interference cancellation process 
is determined by the decoding order which we denote by a 
permutation of the user indices as tt = { 7 ri, 7 r 2 ,-- - 
namely, if tt ^ = k, then Xk (or Xt^) is the i-th user signal to 
be decoded. At each step of SIC, the previously decoded user 
signals can be regenerated by using the same channel coding 
and modulation as having been used by the source and then 
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cancelled out from the received signal. Hence when decoding 
the i-th user’s signal the interference from the j-th user’s 
signal Xtt with j < i can be removed. At the fc-th user , 
the SNR associated with the decoding of the i-th user’s signal 
can be given by 

(3) 

777^0!/ + 1 

where = J2f=i+i denotes the sum of the PAFs of the 
users whose signals are decoded later than Xt^- So 11.,^^ will 
fail in decoding a;,ri if 7^'“ < 2’’”^ — 1. Note that the SNR 
expression in © is based on the assumption that Xt^^, j < i 
have been successfully decoded at the fc-th user. We have the 
outage event of Xt^ at C/,rfc under decoding order tt as follows 

U 7^; <2’'-^ -lLfce/C,z<fc. (4) 

The outage event of the fc-th user is simply the outage event 
of Xt,^ at user i.e., OJJ;. 

Here, our purpose is to find the optimal power allocation and 
decoding order selection to balance the outage probabilities of 
the users under certain fairness constraint, which is formulated 
as the following minimum weighted success probability max¬ 
imization (MinWSP-Max) problem 


max min (1 — Pr { 01 “^)) 

a, 71- kGJC^ 

(5a) 

s.t. TT G n, 

(5b) 

0 < ctTTfc < 1, for fc G /C, TT G n, 

(5c) 

a,rfc < 1, for TT G n. 

(5d) 

ketc 



where is the weighting factor of the fc-th user denoting 
the relative importance of its success probability and 11 is the 
set of all candidate decoding orders. 

When user grouping is considered, the users are divided into 
multiple groups and only one group of users are scheduled in 
each channel block. In other words, the different user groups 
are orthogonally multiplexed. In each group, the transmission 
from the source to the users follows the same way as that in an 
elementary downlink NOMA system. In fact, each user group 
can be treated as an elementary downlink NOMA system with 
fewer users, for which a similar problem as stated in Q can 
be formulated. In order to balance the outage performance of 
the users in different groups, inter-group power and channel 
resource allocation should be considered. The details on these 
problems are deferred to Section HVl 

HI. Solution of the MinWSP-Max Problem 

To solve the MinWSP-Max problem, a main obstacle is the 
decoding order selection which is nonlinear and of combi¬ 
natorial nature. Since the number of the candidate decoding 
orders is K\, which can be large for even moderate values of 
K, it is difficult to acquire a knowledge on the relationship 
between the decoding order and the outage probabilities of 
the users. To circumvent this challenging issue, we first study 
the optimal power allocation problem under a certain assumed 
decoding order. The recognitions acquired on the optimal 


power allocation under a given decoding order make it more 
easier to grasp how does the decoding order affect the fairness 
among the users. Then, by induction, we derive the optimal 
decoding order from the outage balancing perspective. The 
optimal PAFs and balanced outage probabilities of the users 
are obtained in closed-form by investigating the optimal power 
allocation under the optimal decoding order. 


A. Optimal power allocation under a certain decoding order 

Consider any decoding order tt = {tti, 7r2, • • • , ttk}- From 
©, to decode x-^^^ successfully at C/^^, i < k and k € 1C, a 
necessary condition is 


771-fcQ!7i-^ 

1 I ^ ^ 

1 +77rfca7 


( 6 ) 


which can be rephrased as 


7-. [a., - - 1) a^] > 2^-i - 1. (7) 

Obviously, when — (2’^”’ — 1) a^' < 0, the inequality in 
(|7]i can never be satisfied, namely, Xt^ can never be decoded 
irrespective of the channel SNR. Hence, the PAFs should be 
selected such that 


> 0, fc e /C, (8) 


or else the users will always be in outage, which is unwanted. 

In the following, we focus on the decoding of and Xtts 
at an arbitrary user and assume that the PAFs of the other 
signals are given (thus, is known). From (IHi, we have 


={7^„ <7th^} (9) 

= {77r„ < max , 7t?)} - (10) 


where, for convenience, we use the notation 


7th = 




2^-k - 1 


fc e /C 


( 11 ) 


which denotes the channel SNR required to decode given 
that XjTi, i < k have already been decoded and removed from 
the received signal. 

To proceed, we introduce the following proposition on the 
optimal PAFs. 

Proposition 1: At the optimal power allocation, the con¬ 
straint in ( fSdl i is satisfied with equality, i.e.. 


keic 


Proof: When ^ 1’ decrease k G 

K. by scaling up all the PAFs by a factor of e = 1/ 
which gives rise to a lower outage probability for all the users. 
Hence, the optimal PAFs should satisfy (fT2l) . ■ 

From (fT2l i. we know that = 1 — — ctTra ■ Hence, we 

concentrate on the selection of in the following. According 
to = ctjrj + and the constraints in (fT^ and ® (for 
k = 1 and 2), we obtain the feasible range of ctiTj as follows 


a2 e ((2"-= - I)a7^ I/2^-i - . (13) 


It is not difficult to see that as approaches its lower 
bound, turns to be infinite while is a finite positive 
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value, and as approaches its upper bound, 7 ^^^ turns to 
be infinite while 7 ^^^ is a finite positive value. So there exists 
an e (( 2’'-2 - such that 

when a ^2 = Then, by the monotonicity of 7 |^^ and 7 ^^^ 

with respect to we can rephrase in (fTOl i as follows 


O 


TTn, 

7r2 


{l^r. < 7*'} : «2 e (a^, 1/2’'’^! 

{7.„<7r}, a2 G ((2"- -!)«?,«;]. 


(14) 


performance analysis while not only the fairness emphasized 
case as considered in our work. 


B. Optimal decoding order from the outage balancing per¬ 
spective 

Consider the decoding order tt = {tti, 7 r 2 , • • • , and as¬ 
sume that the optimal power allocation conditions in ([ 8 ]), (fT2l) . 
and (flST l are satisfied. Then, the weighted success probability 
of each user U^r^. can be given as follows 

(l-Pr{0;^^})“-'= = [l-Pr| IJ 7 .. <7^^ 

V [fex;.j<fc 
= ( 1 -Pr{ 7 ., < 7 *''}^'= 

= expr-—-Vfce/c (16) 

\ r.n-fc/tC-n-fc / 

where in the second step we have used (fTSl l in Theorem [ 1 ] 
Then, we have the following theorem on the optimal decoding 
order. 

Theorem 2: For the MinWSP-Max problem in Q, there 
exists an optimal decoding order tt ° = {TrJ.Tr",--- , 7 r^}, for 
which the following conditions are satisfied 



Fig.millustrates the variation of the values of 7 ^^ and with 
respect to From Fig.[T] for any a in 1/2’’’"i — 

we can find an a in (( 2’’’"2 — that satisfies 

7 (h''|a„ 2 =a = namely, the required value of 7 ^„ 

to recover is the same when aTr 2 = ci as when = d. 
However, the decoding of becomes easier by setting 
to be d instead of d, since decreases with Thus, 
for the decoding of and Xt^^, it is optimal to select 
in (( 2’'”2 — l)aj^,a* 2 ], which is equivalent to the constraint 
of 7th< 7th^- What is more, to make the decoding of x,ri 
and more easier is beneficial for the decoding of the 
following signals (i.e., x^r^, i > 2). Hence, the optimality of 
the constraint holds true for the whole system. 

The above conclusion can be extended to the selection of all 
the PAFs using an inductive method, for which the following 
theorem is given. 

Theorem 1: To achieve the optimal outage performance for 
a given decoding order tt = ( 7 ri, 7 r 2 ,-- - ■,'Xk), the PAFs 
should satisfy the constraints in ® and (fTST i and 

7r<7r<---<7tr- (15) 

Proof: See Appendix for the proof. ■ 

Theorem [T] implies that with a pre-assumed decoding order, 
the PAFs should be selected such that each user signal is more 
easier to be decoded than the user signals that are decoded 
later than it. From another perspective, if a user’s signal can 
be decoded, then the signals of the users with a prior decoding 
order should be decodable too at the same channel SNR. 
Otherwise, a different decoding order from the one assumed 
could be used to achieve a better outage performance. 

The constraints in ([ 8 ]), (fT 2 l i. and ( fTSl ) form the conditions 
on the optimal power allocation in the elementary downlink 
NOMA system, which can be applied to more general outage 


r7r° r^o 

—- < —^ < 

^ 77 ° W^c 


< 




(17) 


Proof: We consider an arbitrary decoding order tt = 
{ 7 ri, 7 r 2 , • • • jTTif} and prove Theorem |2] by showing that for 
any two adjacent users and 1 < m < AT, if 

> ^TTm+i’ ^hen by exchanging the decoding 
orders of Xt^^ and the minimum of the weighted 

success probabilities of these two users can either be increased 
or keeped unchanged, while not affecting the weighted success 
probabilities of the other users. If the above statement is true, 
then by iteratively optimizing the decoding order of any two 
adjacent users, we can achieve an optimal decoding order that 
is the same as the one in Theorem |2] See Appendix iBl for the 
detailed proof. ■ 

Theorem |2] implies that to balance the outage probabilities 
of the users, it is optimal to assign a higher priority in the 
decoding sequence for the signal of the user with a smaller 
weighted average channel gain, namely, the optimal decoding 
order depends on the ordering of the weighted average channel 
gains of the users, where the weighting factor imposed on the 
average channel gain is the reciprocal of that imposed on the 
corresponding success probability. For instance, \f Ti/wi < 
Fj / Wj , then Xi should be decoded prior to Xj. 

The conclusion of Theorem |2] is intuitive from the weighted 
success probability expression in (fTfil l with the monotonicity 
of the exponential function and the optimal constraints on 7 ^^'“, 
k G 1C in (flSl l. 


C. Optimal PAFs and the Optimized Outage Probabilities 

Suppose that the optimal decoding order as stated in Theo- 
rem|2]is adopted. Then, according to Theorem[T] the MinWSP- 
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Max problem in (|5]) can be reformulated as 


max min exp (--) 

a fee/c j 

(18a) 

s.t. ^ = 1 , 

(18b) 

keK 


0 < a 7 r° < l,for k G 1C, 

(18c) 

> (2''"k — a^'=,for k G 1C, 

(18d) 

7th < 7th < • • • < 7th • 

(18e) 


Since exp(—a:) is a monotone decreasing function of x, the 
objective function in (fTSh ') can be replaced by 
In addition, the constraints in (fTSb-e) can be combined to be 
the following compacted form 

0<7.f <7^ (19) 

Hence, the optimization problem in ( fTSl l is equivalent to 


min 

a 

7th 

max 

keJC F^o/w^c 

s.t. 

E ^ 


fce/c 


constraints in (HSll. 


( 20 a) 

( 20 b) 

( 20 c) 


From the preconditions that Ftt" /wt^o < /wt^o < ■ ■ ■ < 
and the constraints in (I20cb . it is possible to select 
the PAFs such that all k € K. are equal. Inspired 

by this, we have the following proposition. 

Proposition 2: At the optimal solution of ( l20t . the follow¬ 
ing constraints are satisfied 


7th 


7,h^ 


7th'’ 




r,r° /Wti-J. 


= A. 


( 21 ) 


where A is an auxiliary variable. 

Proof: See Appendix |0 for the proof. ■ 

From Proposition |2] the users have equal weighted success 
probability at the optimal solution of the MinWSP-Max prob¬ 
lem, which is exp(—A). The resulted outage probabilities of 
the users are PYf = 1 — exp(A/w 7 r°), k G 1C. As can be 
expected, the user with a larger w will suffer a smaller outage 
probability. The optimal power allocation is given in terms of 
the PAFs in the following theorem. 

Theorem 3: For the elementary downlink NOMA system, 
to achieve the optimal balanced outage performance, the PAFs 
of the users should be selected according to (l 22 l i: on top of 
the next page, where 


A = 


- 1 


K 

E 


2 — 1 r^o 


(23) 


r7rf/t«7rf 

Proof: First, expand and rephrase the equations in (l2Tl i 
as follows 


1 


K 


I P O 


+ E I - 1) , fc € /C. 

j=k+l 

( 24 ) 


Then, from the fact that A = expression of 

A in (l23b can be obtained through some manipulations. The 
optimal PAFs can be obtained by solving the equations in (|2TI) 
successively from large to small with respective to the value 
of k. ■ 

Recall that r.^.^ = P/AoIE(|/i^o p), from the expression of 
a,ro, k G JC, the optimal PAFs are determined by the targeted 
data rates, average channel gains, and weighting factors of the 
users, while not affect by the transmit SNR, P/Nq. Moreover, 
to increase P/Nq only decreases the outage probabilities of 
all the users by the same amount in the logarithmic scale and 
does not affect the relative outage performance of the users. 

IV. NOMA WITH User Grouping 

In the elementary downlink NOMA system, all active users 
are scheduled simultaneously in each channel block. However, 
this may be impractical, because the complexity of SIC scales 
at least linearly with the number of the users that are involved 
in a transmission m, m, M- In this section, we consider 
the case with user grouping, in which the number of the users 
involved in a transmission is much smaller than K, and hence 
a better complexity and performance tradeoff can be achieved. 
The associated problems of power and resource allocation 
among different user groups are investigated from the outage 
balancing perspective. 


A. User Grouping 

To realize NOMA with user grouping, the first problem is 
how to group the users. We focus on the case when all the 
user groups have the same number of users, denoted by L. 
This is considered because it is simple to be realized and (as 
will be shown) can achieve a good enough performance. Use 
0 = {1,2, ■ ■ ■ , G} to denote the set of the user groups, where 
G = K/L. Though it is assumed that AT is a multiple of L, 
the following analysis can be applied to all possible values of 
K. Specifically, when K is not a multiple of L, we can add 
GL—K virtual users with extremely good channel qualities (or 
extremely small weighting factor) to the user set. The system 
with virtual users can well approximate the original one, since 
the power required by the virtual users is nearly zero, which 
has trivial effect on the performance of the other users. 

Two grouping algorithms will be considered, random group¬ 
ing and optimal grouping. At each time, the random grouping 
algorithm randomly selects L users from the ungrouped user 
set and forms them into a new group, which is repeated until 
there are no users remained. The computational complexity of 
random grouping is 0{K). The optimal grouping algorithm is 
realized by recursive search of all candidate growing modes. 
The number of the grouping modes is K\/ ((L!)t which 
increases exponentially with K. 


B. Inter-Group Power and Resource Allocation 

For downlink NOMA with user grouping, each user group 
can be treated as an elementary downlink NOMA system with 
L users. Similarly as in Section UlI-CI there should be an “A” 
for each user group, which we denote by Ag for g G Q. Use 












2""i+i - 1 
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K 


^K+i/'^K+i 

■\T-o —1 




,-n TTro/m 

fc + l/ ^'fc + 1 j^k+2 3' 


J-l 

fc+l’ 


2 - 1 


= 


2 "K " \ 2 - 1 


2'"^K - 1 


r^o Iwtto a 

^K' 


, for fc < - 2, fc e /C, 

, for k = K — 1, 

, for k = K. 


( 22 ) 


Pg and tg respectively to denote the proportion of the power 
and the channel resource allocated to group g. Then according 
to (|2^ . Ag can be given as follows 


2’'si/*3 _ 1 

j\ — _ _ 1 . 

= fgitg)IPg 


^91 1^9 — X 


L 

- 

1=2 


?5Zj=l 


(25) 


where gi is the index of the Z-th user in group g after sorting the 
L users according to Theorem |2] i.e., Tg^/uigj^ < Tg^/wg^ < 


■■ < Tg^Mfc and 


fgitg) = tg 


2 '’bi/‘s _ 1 
^gi/'^gi 


A 2’'9i/‘9 - 
\ —-^- 2^1 = 1 

1=2 


/^9 


^gi/'^ 


(26) 

The corresponding optimal PAFs of the L users in group g 
can be obtained according to (l22l i. which are not shown here 
for save of space. 

Note that the L users in each group g, g & G will have 
the same weighted success probability which monotonously 
decreases with Ag. Hence, the problem of inter-group power 
and resource allocation to maximize the minimum weighted 
success probability of the K users is simply to minimize the 
maximum Ag, g € G, which can be formulated as follows 


min max Ag 

(27a) 

Pg,tg,g£G g&G 


S.t. Pg>0,gGG, 

(27b) 

ypg < 

(27c) 

g&G 


tg > 0,9 G G, 

(27d) 

E^9<1- 

(27e) 

g&G 



For the above optimization problem, we have the following 
propositions 

Proposition 3: At the optimal solution of dZTl i. the follow¬ 
ing conditions must be satished 

Ag = AQ,g G G- (28) 

where Aq is an auxiliary variable. 

Proof: We prove this by contradiction. Suppose that the 
conditions in ( |28] | are not satished at the optimal solution and 
denote by ^ = argmaxggp Ag. It is straightforward from the 
monotonicity of Ag with respect to pg that by decreasing all 
Pg, 9 G 0 \ iy by an appropriate amount 5 and increasing 
Pg, 5 G by (G/G — 1)(5, the value of maxggg Ag can be 


decreased, where G is the cardinality of G- The proposition is 
proved. ■ 

Proposition 4: At the optimal solution of dZTl i. the con¬ 
straints in (I27bl i and (I27db are satished with inequality and 
the constraint in (I27cl i with equality. 

Proof: First, consider the constraints in (I27bl) and (I27dl) . 
if any of them are satished with equality, the corresponding 
group of users will have a zero weighted success probability, 
which violates the outage balancing criterion. Then, if ( I27cb 
is satished with inequality, by scaling up all pg, p G (y by 
a factor of ^IJ^geoPg’ ^ 9 ’ 9 ^ G "'ill be decreased, 
as Ag, 9 G G motonously increases with pg when pg > 0, 
which violates the minimization criterion. Hence, to achieve 
the optimal solution of (l27b . the statements in Proposition |4] 
should be satished. The motonicity of Ag with respect to pg 
is straightforward from (l25b . ■ 

Proposition [3 implies that the K users will have the same 
weighted success probability at the optimal inter-group power 
allocation, no matter how the channel resources are allocated. 
From proposition @1 the constraint in (I27cb is satished with 
equality at the optimal inter-group power allocation. Then, 
according to (l25l l and (l28l l. we have 


Ao — Ao — ^~~^Pg^g — fg^fg)- (29) 

geG geG g&G 

The corresponding optimal pg, 9 G G can be obtained as 

fgjtg) 


Pg = 




(30) 


Now, we consider the optimization of the resource allocation 
parameters, which is further formulated as follows 


min Aq (31a) 

s.t. ^^tg < 1, (31b) 

geG 

tg > 0,9 gG- ( 31 c) 

where the constraints in ( I31cb are based on Proposition |4] To 
solve problem ( l3Tb . we have the following propositions 
Proposition 5: For any p G (y, in the region of {tg\tg > 
0}, the hrst derivative of fg{tg), given in (l32T i. is a strictly 
monotone increasing function of tg, and fg{tg) is a strictly 
convex function of tg. 


L-l 


rg(tg)=-^-yA^- 


^ 9L 


1^1 

In: 


^91+1 


^91 ^9! + l 

^ -Xg) ^9 &G 


cl In xi 


- 4 ) 


( 32 ) 
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where for convenience, we use the notation ig = 2^^=^ . 

Proof: See Appendix |D] for the proof. ■ 

Proposition 6: Problem OTl i is strictly convex. 

Proof: Base on the convexity of fg{tg) (see Proposition 
0 ), it is obvious that Aq is a strictly convex function of tg,g G 
Q in the region defined by (I31cl) . as Aq is a linear combination 
of fg{tg),g G Q. Moveover, the constraints in (I31bl i and (I31cb 
are linear, i.e., convex. Hence, problem OTT i is convex. ■ 
We resort to the method of Lagrange multiplier to solve 
problem OTI) . The Lagrange function is given by 


■^(^1 j ■ ■ ■ AGt ^ ^ fgi^g) I ^ ^ ^g ^ 


(33) 


see 


vsee 


where A is the Lagrange multiplier for constraint (I31bl i. and 
the Karush-Kuhn-Tucker (KKT) optimality conditions are as 
follows 

- = nt 

dtg •'9^ 


see 


ksee 


■A = 0, 

9 ^G, 

(34a) 

tg<^ 


(34b) 

)=o. 


(34c) 

A > 0, 


(34d) 

V 

o 

9 &G- 

(34e) 


Due to complicated expression of fg{tg), it is hard to solve 
the optimality equations in (l34t analytically. Here, we provide 
a simple iterative algorithm to obtained the optimal fg’s, for 
which the following proposition is essential. 

Proposition 7: At the optimal solution of OTI) . the con¬ 
straint in (I31bl) is satisfied with equality. 

Proof: Recall that problem (EB is convex, hence the 
conditions in OTl i are necessary. From 02l i. fg{tg) < 0 always 
holds true for any g G Q in the feasible range defined by (I31cl i. 
Hence, by (I34ab . A has to be strictly positive. Further, (I34cb 
implies that constraint (I34bb or (I31bb must be satisfied with 
equality. ■ 

Since for any g G G, fg{tg) < 0 and fg{tg) strictly increases 
with tg in the feasible range (see Proposition |5]), there exists a 
unique and strict positive solution for all equations in (I34ab for 
any given A > 0. With this observation, the problem of finding 
the optimal fg’s can be alternatively solved by finding the A 
which yields fg’s that satisfy (I31bb with equality. Though it is 
hard to find fg’s adding up exactly to one, we can iteratively 
bound A to mitigate the gap of X^gep ^g upper bound 

to an acceptable value. To realize this, we have the following 
algorithm 

Step 1; Set Xl = 0, Xh = — maxggg /g(l), 0 = 0; 

Step 2; Set A = {Xl + Ac/)/2 and solve (I34ab to obtain tg, 
g &G; 

Step 3; If Yhg^g tg <1, 9 = 0, else 0=1; 

Step 4: If Xh - Xl < Co, t* = tg for g G G; 

else, set A^, = {1 — 9)Xl + 0X and Xh = OXh + (1 — 
0)A, and go to Step 2. 


TABLE I 

Randomly generated r, w, and d for the results in Fig. [2] 



r (bits/s/Hz) 

w 

d (meter) 

G! 

(0.57, 0.04, 1.39) 

(0.39, 0.30, 0.31) 

(0.453, 0.788, 0.417) 

G2 

(0.91, 0.35, 0.74) 

(0.11, 0.29, 0.60) 

(0.535, 0.981, 0.480) 

G3 

(0.47, 0.74, 0.79) 

(0.25, 0.35, 0.40) 

(0.904, 0.842, 0.208) 

G4 

(0.65, 0.23, 1.12) 

(0.45, 0.46, 0.09) 

(0.636, 0.550, 0.870) 

G5 

(0.73, 0.22, 1.05) 

(0.32, 0.26, 0.42) 

(0.951, 0.531, 0.784) 


Note that from Step 2 to 4, many other line search algorithms 
can be used. In step 2, the equations in ( I34ab can be solved 
using the Newton’s method or many others ll27ll . The compu¬ 
tational complexity of the algorithm depends on the precision 
required. Use Ci to denote the precision required when solving 
(I34ab , then the total complexity is GlogCologei. If parallel 
computing is available, the complexity will be logCologei. 

The optimal Aq and pg’s can be obtained by substituting 
tg’s into ( l29b and (l30b . respectively. 

V. SIMULATION RESULTS AND DISCUSSION 

In this section, we present simulation results to verify our 
analysis and show the potential performance gain of NOMA 
over OMA. The tradeoff between the complexity and perfor¬ 
mance of NOMA is also discussed. We consider a network 
model where the users are uniformly distributed over a circular 
disk (with its radius being normalized to be 1 meter) centered 
at the source. We assume that Hk = dk~'^ with dk being 
the distance from the source to user Uk and g the pathloss 
attenuation factor which is set to be 3.75. The other parame¬ 
ters, including the targeted data rates and weighting factors of 
the users, will be specified for each of the following results. 
For convenience, we use the notations r = {ri, r 2 , • • ■ , rjf}, 
w = {wi,W 2 , ■ ■ ■ ,wk}, and d = {di, d 2 , • ■ • , As 
a reference scheme, we use the conventional TDMA which 
refers to the orthogonal allocation of the channel resources and 
is essentially equivalent to any OMA scheme uni Sec. 6.1.3]; 
both the cases with only power allocation (PA) and with joint 
power and resource allocation (PARA) are considered. The 
optimal resource allocation in a TDMA system can be solved 
by using a methodology similar to that for inter-group resource 
allocation in NOMA (see Section HVl i. 

A. Verification of the optimal power allocation and decoding 
order selection 

To substantiate our analytical results, we provide the solu¬ 
tion of the MinWSP-Max problem by using exhaustive search 
of the optimal decoding order and PAFs in Fig. |2] where we 
take the case of AT = 3 as an example, and randomly generate 
five group (G1 to G5) of simulation parameters, which are 
given in Table U As can be observed from Fig.|2] the results by 
using exhaustive search well consolidate our analytical results, 
in terms of both the outage probabilities and the PAFs of the 
users. Also, as analyzed, the user with a larger weighting factor 
will have a smaller outage probability. 

B. User fairness enhancement by NOMA 

To evaluate how does NOMA performs on the average, we 
simulate the elementary downlink NOMA system for 5000 


















Fig. 2. Outage probabilities and PAFs of the users using both the analytical 
results and exhaustive search under five group of randomly generated sim¬ 
ulation parameters as given in Table Q] The case of X = 3 is considered, 
the sum rate is set to be = 2 bits/s/Hz, and the transmit SNR is 

P/No = 10 dB. 



Fig. 4. The average outage probability of user 1 with respect to r-^ when 
P/No = 15 dB and K = 10. 




Fig. 5. The transmit SNR required by NOMA and TDMA to achieve an 
Fig. 3. The average outage probability of the user 1 with respect to K when probability of 0.1 for the users when K = 10. 

P/No = 15 dB and = 3 bits/s/Hz. 


times to cover a large number of scenarios; then, for each 
user, the average of its outage probability over all simulations 
(average outage probability in short) is used as a performance 
measure. At each time, the simulation parameters r, w, and d 
are generated independently. The data rates follow a uniform 
distribution under the constraint that ''’fc = where rs 

is the sum of the targeted data rates of the users, which can 
also be referred to as the targeted system spectral efficiency. 
The weighting factors follow a uniform distribution over [0,1]. 
By symmetry of the users, they will have very close average 
performance. Hence, we only need to focus on one of the K 
users in the plots. Here, we choose user 1. 

Fig. [3 shows the average outage probability of user 1 with 
respect to K. It can be seen that the performance gain of 
NOMA over TDMA increases with K when K is small, and 
turns to be steady as K becomes large. This can be explained 
by the fact that SC utilizes the diversity in the channel qualities 
of the users as a new freedom for potential performance gain 


m . Generally, more users implies more rich channel diversity, 
while as K becomes large, the increment will be less evident. 

Fig. a shows the average outage probability of user 1 at 
different sum rates. Fig.|5]shows the transmit SNR required by 
different schemes at different sum rates and at an average out¬ 
age probability requirement of 0.1 for the users. It is observed 
from both the outage performance and power consumption 
perspective that NOMA always performs better than TDMA, 
no matter when PA or PARA is adopted by TDMA. Also, we 
see that the advantage of NOMA becomes increasingly more 
evident as the data rate increases. This is because NOMA 
exploits the channel resources more efficiently than OMA, and 
hence is more benehcial at higher data rate. Note that when 
the data rate is rather high, the advantage of NOMA will be 
less evident, because it becomes difficult for all the schemes 
to support the communication, and both NOMA and OMA 
will have a bad outage performance. The results are omitted 
here for save of space. 
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Fig. 6. Impact of user grouping algorithm and group size on the performance of NOMA. P/Nq = 15 dB and K = S. 


1) Discussion: From Fig. |4] and |5] the performance gain 
of NOMA over TDMA is limited when resource allocation 
among the users is available for TDMA. However, the effect 
of resource allocation in practical systems will be much worse 
than it theoretically does. Theoretically, the channel resources 
can be arbitrarily divided among the users, i.e., continuous 
resource allocation. Such an ideal condition is not available 
in practice. Generally, according to the operating mode of 
existing communication systems, we can only allocate integer 
resource blocks to each user, i.e., discrete resource allocation. 
In this sense, resource allocation will result in additional delay, 
for which the range of resource allocation will be restricted. 
Even if we divide the original resource block into smaller 
parts, it is the most possible to divide it into equal parts. 

For better understanding, we also show TDMA with discrete 
time allocation in Fig. 0] and |5] where we use T to denote 
the number of sub-timeslots that each original timeslot can 
be divided into. Obviously, T reflects the ability in resource 
allocation of the practical system. As discussed above, discrete 
resource allocation performs much worse than the theoretically 
optimal one. 

C. Performance and complexity tradeoff 

In this subsection, the case with user grouping is evaluated. 
The focus is on the complexity and performance tradeoff of 
NOMA. For notational convenience, we denote by NOMA the 
case without user grouping, NOMA-O the case with optimal 
user grouping, and NOMA-R the case with random user 
grouping. In addition, we use PA to denote that only inter¬ 
group power allocation is adopted by NOMA and use PARA 
to denote that both inter-group power and resource allocation 


are adopted by NOMA. The configuration of the parameters 
r, w, and d is the same as in the previous subsection. Also, 
we focus on user 1 in the plots. 

Fig. |6] gives the average outage probability of user 1 with 
respect to rs under different user grouping algorithms and 
group sizes. From Fig. |6] it is observed that both NOMA-R 
and NOMA-O outperform the corresponding TDMA scheme 
(either with PA or with PARA), irrespective of the group 
size L. Though NOMA-R performs worse than NOMA-O, 
the performance gap is small, especially when L is ralatively 
large and/or resource allocation is adopted. Also, we see that 
NOMA-R with L = 2 can achieve a performance very close to 
that of NOMA when PARA is adopted. Even when resource 
allocation is unavailable, NOMA-R with L between 2 and 4 
can reap a large portion of the performance gain of NOMA 
(over TDMA with PA). These obvervations imply that NOMA 
with rand user grouping and a small group size can achieve 
a good enough performance. We know that the complexity of 
inter-group resource allocation in NOMA is the same to that 
of time allocation in TDMA, hence the additional complexity 
of user grouping mainly lies in the user grouping algorithm. 
Since the complexity of random user grouping is negligable 
compared to that of SIC, it is attractive and also effective to 
improve the compexity and performance tradeoff of NOMA 
by user grouping. 

In the previous subsection, we have discussed that due 
to the physical limitations in practical communication sys¬ 
tems, resource allocation can only be performed in a discrete 
manner, which may degrade the effectiveness of resource 
allocation. Pig. |7] illustrates the impact of such restriction 
on the effectiveness of inter-group resource allocation in a 
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Fig. 7. Impact of discrete resource allocation on the performance of NOMA with user grouping. P/Nq = 15 dB and K = 12. 


downlink NOMA system. Interestingly, the case with discrete 
resource allocation and T = 3 achieves a performance very 
close to that with the optimal continuous one. This is because 
in NOMA, the performance of the users in the same group 
have already been balanced by nonorthogonal transmission 
with fairness power allocation. As a result, the imbalance in 
the performance of different user groups becomes less severe, 
and hence the demand for resource allocation (to enhance user 
fairness) is much weaker than in a TDMA system. 


VI. CONCLUSION 

The outage balancing problem was investigated to achieve 
the optimal fairness outage performance in a downlink NOMA 
system when only statistical CSI is available at the transmitter. 
Both the problems of power allocation and decoding order se¬ 
lection were thoroughly studied and solved analytically. It was 
proved that the optimal decoding order is determined by the 
ordering of the weighted average channel gains of the users. 
The simulation results for both the cases with and without user 
grouping illustrated that NOMA performs much better than 
OMA in terms of the fairness outage performance, especially 
when the practical discrete resource allocation restriction is 
taken into account. For the case with user grouping, we solved 
the problems of inter-group power and resource allocation. It 
was substantiated by simulations that user grouping serves as 
an effective method to reduce the implementation complexity 
of NOMA due to SIC. The complexity issue of user grouping 
can be well circumvented by using a random user grouping 
algorithm, which has negligible complexity when compared 
with that of SIC. It was shown that NOMA with random user 


grouping and a small group size (2 to 4 users in each group) 
can reap most of the performance gain of NOMA. 


Appendix A 
Proof of Theorem[T| 

It has alreadv been nroved in Section IIII-Althat < 7t?, 
we only need to further prove that 7th'" < 7*'"^' for 2 < m < 
K under the assumption that 

<ll^ (35) 


Consider the decoding of any two adjacent signals and 
2 < TO < Ff at an arbitrary user 2 < n < K 
under the assumption that the constraints in ([8]l, (fTSI) . and dTSl) 
are satisfied and the PAFs of all the users other than and 
are given. Then, the outage events in decoding 
and a;,rm+i at C/ 7 r„ can be given by 


Ol-_ y 7^; < 2^-. - 1 


and O. 


TTm + l 



U 77r„ < 7t? ? 


yi^K.j<rn J 


= { 7 ^, 

. <7tr} 

(36) 

= [It,, 



= { 7 ^, 

. <niax(7,;™+\7^^-)}, 

(37) 


respectively. It is worth mentioning that the second and third 
steps in (l3^ are based on the assumptions stated in ® and 
drSl) . respectively. 
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Based on = 1’ of/’" = c^u-m+i + and the 

constraints in (l 8 ]l (for k = m and m + 1 ), the feasible range 
of is obtained as 

a^m+i e 1^(2’'’''"+! - 1 )q;^”+ 

(38) 

Following the same lines as the discussion on the selection 
of ofTTa in section IIII-AI it can be prove that for the decod¬ 
ing of a^TTm+i, and the following signals (i.e., 
a^TTm+g, • • •, a^TTfc), ths optimal and a,rm+i should satisfy 
7 *"“ ^ 7 *™^'^- Note that the selection of and 
does not affect the decoding of the signals prior to (i.e., 
Xtti j a;,r 2 ) • ■ ■ 5 it can be concluded that the constraint 

7 ^™+! ^ jg optimal for the whole system. Since and 
are arbitrarily assumed, the optimality of the constraints 
in (flSl l can be proved. 


-| v-^m— 1 

1 ^ ~ z^k^l 


*‘'^k TTm + l 

- ar 


Appendix B 
Proof of Theorem[2] 


Assume optimal power allocation among the users accord¬ 
ing to Theorem[T] Then from < 71 ^”“^^ in (flSl l. we obtain 
an upper bound on as follows 


D^TTm + l — 


2 ’'-m + l _ 1 

+ l — 1 


(1-a), 


(39) 


where a = Efceic,fe#m.m+i denotes the sum of the PAFs 
of all the users except and First, we will show 

that for any feasible values of and we can find 

new PAFs and for and respectively, 

such that the following constraints are satisfied 


“t“ /^7r,n + l — ^ tT, 

7 r^^ < 7 r < 7 r■"^ 

where in (HTt . we use the following notations 

2 ’'’'m + l _ 1 

7th'"'"' = 


2 ’’-™ - 1 


7th = 


In fact, the first inequality in (HTI) can be satisfied if 

2 ’'-- - 1 


> 11 - 


(40) 

(41) 

(42) 

(43) 


(1 - a), (44) 


2 ’--r„ 2 '’'-+l - 

which is definitely feasible, and the second inequality in (HTt 
is equivalent to 


/37t„+i < 1 - a - 


2 ’'-™ - 1 
2 ’--m+i _ 1 


^TTm + l’ 


(45) 


which is feasible without violating the first one considering the 
precondition in (|39] |. The inequalities in dTO imply that if we 
use the new PAFs for XT^m and Xtt^ should be decoded 

after while the decoding orders of the signals 

Xnm+s’ ■■■’ "'dl not be affected. The equality in (l4Ql i 
implies that the PAFs of the other users are not changed. So if 
we keep the decoding orders of the signals prior to x-j^m (i e., 
■ ■ ■ ’ ^ttk) unchanged, the outage probabilities 


of all the users other than and will not be affected 

by the selection of the new PAFs. 

Now suppose that we have exchanged the decoding orders 
of x^^ and X 7 r,„+i by using the new PAFs and 
that satisfy (l44l i and (l45T l, and the decoding orders of the other 
user signals are unchanged. Use and OZZXi 1° denote 
the outage events of user t/,rm and , respectively, under 

the new decoding order. Then, we have 


= {Ttt^+i < max {jXir+^, 7th'"”')} . (46) 

o:: = { 7 ^,„ < 7 r}U^-:+i 

= < iXC } U {3'^- < (^*'"^' > 7th'"”') } 

= < max > 7th'"”')} • (47) 

Note that in the second step of (l47l i. we have expanded O^Z+i 
which can be similarly treated as OZZX\ and in the 

third step, we have used the inequalities in (lUT i. From (l46l l 
and (l47l i. we obtain the weighted success probabilities of C/tt,,, 
and under the new decoding order as follows 


(i-pr{o:::;}) 


= exp 


(i-p,{ 0 wt}) 


= exp 


max( 7 ,^-, 7 ;;-‘ ') \ 

^TTm J 

(48) 

rTTm + l/'a^TI’m + l j 

(49) 


Under the original decoding order, the weighted success prob¬ 
abilities of C/,rm and UT^m+i are given by 


(l-Pr{0:-})”-=exp 
(l-P,{0-; })”-"*■= exp 


T^Tn \ 

7th \ 

TTm + l 


r-TTm + l/^TTm + l 


(50) 


(51) 

With the facts that 7 ^'"’*"'^ < 7th™ < 7*'"^' 7th'"”' — 

7 th™ ^ 7 *""^'’ proved that when > 

r-7r,n + l /ra-TTm + l ’ 


-in((i-pr{a::})“-,(i-Pr{a-j})“-) 

> min ((1 - Pr {O:-, (l - Pr . 

(52) 

In other words, if r,r^/rt; 7 r„ > r 7 r,„+t/iU 7 r,„+i, the minimum 
of the weighted success probabilities of and UT^m+i after 
exchanging the decoding orders of their signals is greater than 
or equal to that under the original decoding order. 


Appendix C 

Proof of Proposition!!] 

We will prove this by controdiction. Suppose that not all 
7 th"• Wtt" are equal at the optimal solution and denote 
by K,' the set of the indices such that 7 **''/r,r°^ • = 
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maxfcgyc lib ^ *^P 

a,r°,, k' € 1C' hy 6' > 1 and scale down all a^o, k £ IC\IC' hy 


6 = 


^k^K\K' ' 


JT° J2k'GIC' 


.?,(l--5') 

fc' _ 


< 7t?Vr.o 


for 


^ ^^^/ CK^o ^ ^tn / 

all A:' £ 1C can be decreased while 7 ^^*“/F^o • , k € IC\IC 

will increase, which is straightforward from the definition of 
7 th (see (fTTT i). The value of S ensures that the constraint in 
(l20b l is satisfied. In addition, it can proved that if S' is small 
enough, the constraints in (l 20 b ) will not be violated, and the 
maximum of /T^^o ■ Wt^o, k € IC\IC can be smaller than 

the minimum of 7 (^'“'/r 7 r“, • tr'ii-", > k' £ 1C. The proof of the 
second part is stratightforward, hence we focus on the proof 
of the first part in the following. 

Denote by 7 ^^'“ the value of 7 ^^'' after scaling up or down 
the PAFs as stated in the above. Then, according to (fTTT i. the 
expression of 7 j^'“ is given in (|53]i; on top of the next page. 
First, we prove that there exists a <5' > 1 satisfying > 0. 
With the fact that 7 j^^ > 0 , it is directly from (fTTl i and (I53bb 
that 7 j^^ > 0 is always satisfied if 1 G 1C. While if 1 £ 1C\1C, 
7 (^^ > 0 can be satisfied only if we choose a small enough 
S' {> 1), since > 0. Use S'l (> 1) to denote the upper 
bound on the value of S' that satisfies 7 ^^^ > 0. Then, we 
prove that there exists a b' £ ( 1 ,^ 5 .] satisfying 7 jjj'“ < 
for the following four cases. 


• Case 1: k € 1C\ 1C, k + 1 G 1C\1C, and k < K-, 

• Case 2: k G tC, fc + 1 £ 1C, and k < K', 

• Case 3: fc £ 1C, fc + 1 £ /C \ 1C, and k < K', 

• Case 4: k G 1C\ 1C, fc + 1 £ 1C, and k < K. 

It can be proved that is always satisfied for the 

7r° 

first three cases by simply comparing the expressions of 
and with the fact that 7 ^'' < 7 ^*'^F For Case 4, we 

have the following relation of equivalence 



2''-l - 1 


S - S 2 ’''^+! - 1’ 
(54) 


which can be easily obtained from (l53T l. From 7 ^'' < 
we know that 



> 


- 1 


(55) 


Obviously, if the condition in (l55l l is satisfied with inequaltiy, 
we can find a > 1 such that (l54l i is satisfied if S' < S' 2 . 
It can be seen that the condition in (l55l l can only be satisfied 
with inequaltiy in Case 4. Since if not so, 7 ^*“^^ will 

be equal; then from the fact that Ftt "/wtto < Ftt" 7 ^,^“ , 

k ' k fc + 1 ' fc + 1 

1th will be greater than or equal to lib'^^/T^o^^ ■ 

which violates the fact that k G 1C\ 1C and fc + 1 G 1C 
in Case 4. From the above discussions, if we choose S' in 
the range of ( 1 , min ^ 2 )]’ '^he constraint in (l 20 b ) will be 
satisfied. 

It follows that if we choose a suitable S', max^gx: 1^'° ' 

Wt^o can be further decreased, which contradicts the optimality 
assumption. 


Appendix D 

Proof of Proposition[5] 

We first consider the first derivation of fg{tg), which is 
given by 


m,) = 


2^91/tg _ 1 _ 2^91/t9 ln2rgjtg 

Tgi /u’gi 


2^92 1^9 - 1 

^52 7^52 


1^92 7^9 

2^911*9 ln2rg^/tg 


2-91C9 _ 1 _ 2 ''nC 9 \n 2 rgjtg ^yi-_t 


r9,7' 


91 / ^91 


1-1 


2'^9i/^9 — \ i__i r /t 1 / 

-^- 2l^j=i^9jtt9 ln2 > ra-/ta 


2r9j*9-l-2-9j*9 \n2rgJtg ^yL-i 




9L / ^91 


L-1 


P In 2 ^ Tg^ Itg. (56) 

^9Ll'^gL 


After some algebraic manipulations, the expression in ( l56l l can 
be rephrased as in ( |3^ . 

Since a; In a; — x strictly increases with x when a; > 1, 
and for any I G C = {1,2, ••• , L], xi > 1 and xi strictly 
decreases with tg when > 0 . It is then clear that each term 
x^g Inx^g — x^g, l G L K dt stactly monotone decreasing function 
of tg when tg > 0. In addition, it is always true that Wg^ /Tg^ > 
0 for ( £ £ and WgjTg^ > WgjTg^ > ■ ■ > WgjTg^. 
Hence, f'g{tg) is a strictly monotone increasing function of tg 
when tg > 0 . 

With the knowledge that f'g{tg) strictly increases with tg 
when tg > 0 , it is intuitive that the second derivate of fg{tg) 
will be positive when tg > 0 . However, for completeness, we 
still give the second derivate of fg{tg) as follows 


l^l 91 ^ 9l-\-l / 

+ (57) 

^ 9L 

It is obvious that fg{tg) > 0 when tg > 0 , which proves the 
strict convexity of fg{tg) with respect to tg when tg > 0 . 
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